Factorization theorem for high-energy scattering near the endpoint 
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A consistent factorization theorem is presented in the framework of effective field theories. Con- 
ventional factorization suffers from infrared divergences in the soft and collinear parts. We present 
a factorization theorem in which the infrared divergences appear only in the parton distribution 
functions by carefully reorganizing collinear and soft parts. The central idea is extracting the soft 
contributions from the collinear part to avoid double counting. Combining it with the original soft 
part, an infrared-finite kernel is obtained. This factorization procedure can be applied to various 
high-energy scattering processes. 



Theoretical predictions on high-energy scattering pro- 
cesses are based on the factorization theorem, in which 
the cross section is factorized into hard, collinear and soft 
parts. The hard part consists of the partonic cross sec- 
tions with hard momenta. The collinear part describes 
the effects on the energetic particles. For hadron col- 
liders, the collinear part involves the parton distribution 
functions (PDF). The soft part describes the soft gluon 
exchange between different collinear parts. These three 
parts are decoupled and the collinear and soft parts are 
defined as the matrix elements of gauge-invariant opera- 
tors. 

This conventional factorization scheme, as it stands, is 
plagued by infrared (IR) divergences in the collinear and 
soft parts. The hard part, by definition, depends on hard 
momenta and it is IR finite. There may exist ultraviolet 
(UV) divergences, but they can be removed by countert- 
erms. The collinear and soft parts unavoidably encounter 
IR divergences in massless gauge theories like QCD what- 
ever the regulators are. If the dimensional regularization 
is employed to regulate UV and IR divergences, the di- 
vergences are expressed in terms of the poles in euv ( £ ir) 
for UV (IR) divergences. These UV and IR poles appear 
separately and there is no mixing between them. If the 
UV divergence is regulated by the dimensional regular- 
ization, but the IR divergence is regulated by the offshell- 
ness of external particles, there are UV poles, and also 
logarithms of the offshellness for IR divergences. 

Our idea of factorization goes further than the opera- 
tor definition and its radiative corrections of each part in 
order to isolate IR divergences. The radiative corrections 
are reorganized such that the IR divergences reside only 
in the PDFs, as in full QCD, and the remaining quantities 
are IR safe. Furthermore the method of the reorganiza- 
tion is so simple and universal that it can be applied to 



a variety of scattering processes. The goal is to find an 
IR-finite kernel to be convoluted with the hard function 
and the PDF in the factorized expression of the cross 
section. The evolution of the kernel can be consistently 
described using the renormalization group equation with- 
out devising a method to handle IR divergences. We will 
delineate our method of factorization, and explain how 
it describes the factorization theorems in deep-inelastic 
scattering (DIS), and Drell-Yan (DY) processes near the 
endpoint or with transverse momentum. 

We can schematically write the factorized scattering 
cross sections near the endpoint as 



da = 



H(Q, n) ® S(E, n) ® fi(E, p) ® f 2 (E, n), 
H(Q, p) ® J(QVT^, fx) ® S{E, n) ® f(E, p), 

(1) 

in DY and DIS processes respectively. H is the hard 
function, and S is the soft function, which is defined in 
terms of soft Wilson lines [l} . fx, f% and / are the initial 
collinear functions denned in terms of the quark bilinear 
operators at the intermediate scale E ~ Q(l — z). J 
is the jet function in DIS, which describes the collinear 
final-state particles where z is the Bjorken variable. Here 
® means an appropriate convolution, and Q is the large 
scale in the scattering. The factorized form in Eq. ([1]) is 
well established in QCD HQ. 

Since the processes involve multiple scales, it is con- 
venient and transparent to understand the physics by 
employing an effective theory, which is soft-collinear ef- 
fective theory (SCET) [4, 5]. The conventional factoriza- 
tion theorem, Eq. ([1]), can be rephrased exactly by noting 
that SCET factorizes the soft and collinear interactions 
from the outset @,0|- The soft part is expressed in terms 
of the soft Wilson lines as 

S OY (E,fj,) = l-(Q\ixY^Y n 5(l-z+2ido/Q)Y^Y n \0) (2) 



2 



for DY process near the endpoint z — Q 2 /s ~ 1 where s 
is the partonic center-of-mass energy squared. And 
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\tiY^Y n S(l-z+in-d/Q)YiY n \0) (3) 



for DIS near the endpoint z = Q/n ■ p ~ 1 , where z is the 
partonic Bjorken variable. The typical soft scale is given 
by E ~ Q(l — z). Here n and n are lightcone vectors 
satisfying n 2 = n 2 = 0, n ■ n = 2. Y"„ and Y fi are the soft 
Wilson lines 0, @] given as 



y„ = exp 
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where A s is the soft gauge field, and V is the opera- 
tor extracting the momentum from the gauge field. The 
one-loop corrections for the soft parts are given, using 
the dimensional regularization for both UV and IR di- 



vergences, as 
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As can be seen clearly, the soft functions involve both UV 
and IR divergences, therefore the evolution of the soft 
function at this stage does not make sense. Furthermore, 
the term 1/(1 — z)+ definitely comes from the real gluon 
emission. In DY process, it is IR divergent contrary to 
the belief that it originates from UV divergence. 

In employing SCET, it is critical to separate the de- 
grees of freedom into collinear and soft parts carefully. 
Special care should be taken in considering loop correc- 
tions in the collinear part since the loop momentum can 
reach the soft region which should be removed from the 
collinear part. In SCET, the zero-bin subtraction [9j is 
devised to remove the soft modes from the collinear part 
to avoid double counting. The idea is to take each Feyn- 
man diagram for the collinear contributions, and subtract 
the contribution when the loop momentum becomes soft. 
Only after the zero-bin subtraction, the separation of the 
collinear and soft parts becomes complete. We emphasize 
that the zero-bin contribution should have commensurate 
scales with the soft part. This needs some elaboration. 
There is a hierarchy of scales Q 3> E ^> Aqcd- We first 
construct SCETt between Q and E and the soft zero-bin 



contribution is that of the collinear part with momentum 
scaling as QX ~ E. There might be the ultrasoft (usoft) 
contribution of order QA 2 to be subtracted, but it is not 
the right choice. If we naively subtract the usoft modes 
from the collinear part, there is mismatch between the 
soft part in Eqs. ©, ^ and the zero-bin subtraction. In 
SCETn below E, we match each quantity at the bound- 
ary E, and evolve down to some renormalization scale /i. 
Here in SCETn, we also perform loop corrections, but 
the zero-bin contribution here corresponds to the usoft 
contribution of order QX 2 . 

The initial collinear function / in SCETt is defined as 

f{x,n) = (N\x n '-S(Qx-n-P)xn\N), (7) 

where \ n = W^ n is the collinear gauge-invariant 
fermion with the collinear Wilson line W n given by 



Wn = exp 
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(8) 



The S in Eq. © is the regulator to treat rapidity diver- 
gence, if there is any. It should be matched to the PDF 
<j> in SCETn as 



f dz fx \ 

f(x, fi) = J — K ( z , m) ■ 



(9) 



The PDF 4> at the operator level is the same as Eq. ([7]) 
for /, but at the scale /j, -c E. The matching coefficient 
K (z, fi) can be computed by comparing different kinds of 
zero-bin subtractions in SCETt and SCETn. 

Now we explain the distinction between the soft and 
usoft zero-bin contributions. In Eqs. ((2]) and (j3|), the soft 
functions include delta functions. The soft contribution 
means that ido/Q or in ■ d/Q is of order 1 — z, and the 
delta functions play a nontrivial role in loop computa- 
tions. If ido/Q or in ' d/Q is much smaller than 1 — z, it 
can be neglected, and the delta function becomes 6(1 — z). 
The remaining soft Wilson lines cancel to yield 1. The 
usoft contribution corresponds exactly to this case. 

The distinction also applies to the collinear function / 
and the PDF cj>. In Eq. ([7]), there also appears a delta 
function. Now consider the zero-bin contribution, in 
which the loop momenta become small compared to the 
label momenta. Near the endpoint, the leading momen- 
tum almost cancels. Therefore in real gluon emissions, 
the small subleading momentum of order QX should be 
kept in the delta function. The soft zero-bin contribution 
is the soft limit of the collinear loop calculations of order 
QX. For </>, we take the limit where the loop momentum 
approaches QX 2 , and the usoft momentum is neglected 
in the delta function. Therefore the usoft zero-bin con- 
tributions are always proportional to 5(1 — z). 
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At one loop, at the parton level is given by 
/«(*,/*) = (10) 
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where J^ 1 - 1 is the naive one- loop result for Eq. ([7]), and 
/ s q is the soft zero-bin contribution. The PDF 



the parton level is given as 
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where is the usoft zero-bin contribution. Note that 

= and /$ = - /W. The PDF obeys 
the standard DGLAP evolution equation from the UV 
behavior of Eq. (fTTj) . Near the endpoint for DIS and DY 
processes, the initial-state jet function K^- 1 ' at one loop 
is given by 

_ a s C F / 1 1 \ 
7r V euv eiR ' 
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Since the collinear contributions for / and <f> are the same, 
K actually comes from the difference between the soft 
and the usoft zero-bin contributions. Furthermore, since 
the usoft zero-bin contributions vanish, K comes from the 
soft zero-bin subtraction. Note that also includes 

the IR divergence. 

The final-state jet function J is IR finite, and contains 
only the UV divergence. Including the zero-bin subtrac- 
tion it has been computed in [loj ]. The one-loop results 
show explicitly that the collinear and soft parts except 
the final-state jet function include IR divergences, which 
hinders the factorization theorem. Therefore the oper- 
ator definitions of the soft and collinear parts are not 
enough to guarantee the factorization. 

We propose to combine the initial-state jet functions 
and the soft function to obtain an IR finite kernel. In 
SCETu below E, the scattering cross sections in Eq. ([1} 
can be expressed as follows: For DY process, it is 
schematically written as 



da D Y = H(Q, n) <g) S(E, n) ® Ky.{E, fi) 
®K 2 (E,n)® <j>! (/x) <8> 4> 2 (a*) 
= H(Q,fj,) ® W B y{E, At) ® Mp) ® ^(m), 
where the kernel Wjyy is defined as 



(13) 



W DY =S(E,ii)®K 1 (E,n)®K 3 {E,ii). (14) 



For DIS, it is written as 



dams = H(Q, n) ® J{Qy/l - z, fi) ® S(E, At) 

®K(E,n)®ct>{ii) (15) 
= H(Q, At) (8 J{QVl-z, At) ® Wdis (E, fi) ® </>(At), 



where the kernel Wdis is defined as 

Wbis=5(E,/i)(8A'(E,/i). 



(16) 



And the hard factor is given by H(Q,n) — |C(Q,/i)| 2 , 
where C(Q,/j,) is the Wilson coefficient of the back-to- 
back current in SCET. 

These kernels are free from IR divergences, and can be 
computed order by order in perturbation theory. From 
the viewpoint of SCET, the kernels are the matching co- 
efficients between SCETi and SCETn, free of IR diver- 
gences, and should depend only on Q and E, or 1 — z. 
And the renormalization group behavior can be probed 
for the hard part, and the kernel without employing com- 
plicated methods to handle IR divergences. To one loop, 
the kernels are given as 



Wjyy = S(l-z 
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(1 - Z)+ V<| \ 

Wdis = S(l - z) 



Eqs. (|13l) and (|T5|) are our factorization formulae for 
the DY and DIS processes near the endpoint. As in the 
conventional factorization theorems, we first start with 
the factorization of the hard, collinear and soft parts. 
The collinear and soft parts are defined in terms of the 
matrix elements of gauge-invariant operators. Due to the 
presence of the IR divergences in each part, the collinear 
and soft parts are reorganized in such a way that the 
kernel W and the collinear part J in DIS are infrared 
finite. In this procedure, the matching coefficient of / 
and the PDF between SCETi and SCETn is included in 
the kernel W. In the form of Eqs. (JT3J) and ([15]), the 
kernel and the hard part are physically meaningful and 
the evolution of these can be considered unambiguously. 

If a process involves small transverse momentum, the 
relevant collinear function satisfies 



f(x,k ± ,n) 



1 dz 



-K 



.(«,k x ,M)0(^.A*)- (18) 



Combining K±_ with the transverse-momentum- 
dependent soft function, we can also obtain an IR 



4 



finite kernel W± in a similar way Woy is obtained 
near the endpoint [ll|. In this case, there appears an 
additional divergence called the rapidity divergence. 
This rapidity divergence can be regulated by a rapidity 
regulator [12j, and it is shown that the rapidity diver- 
gences cancels in each collinear sector if the regulator is 
introduced only in the collinear Wilson lines 11|. It is 
also true near the endpoint. 

Remarkably the kernel for DIS becomes 5(1 — z) to 
one loop since 

S uis + = °' and the inclusive DIS 
cross section near the endpoint consists of the hard part, 
the jet function and the PDF cf> at the scale fi < E. It 
means that the soft contribution and the initial-state jet 
function cancel. Therefore the PDF satisfies the ordi- 
nary evolution equation, as in full QCD. If we separate 
the soft contribution and the collinear contribution, the 
corresponding collinear part f = K ® <f> does not satisfy 
the DGLAP equation, which was discussed in previous 
literature @, [lfl 0. The fact that Wdis = 5(1 - z) 
is shown to be true to all orders in a s by a simple ar- 
gument. The usoft zero-bin contributions vanish in (f>. 
The soft zero-bin contribution from / is obtained by in- 
tegrating out the momenta of order QX. This is exactly 
the procedure to obtain the eikonal form of the soft Wil- 
son line in the n direction. Therefore the initial-state jet 
function, which is the negative value of the soft zero-bin 
contributions, always cancels the soft function at higher 
orders in a s . In contrast, it is different in DY processes 
since the soft function involves an interaction between n 
and rl-collinear fermions, while the collinear part which 
interacts only within each collinear sector does not pro- 
duce the same soft interaction in the zero-bin limit. 

Now that the explicit results of factorization are pre- 
sented, the explanation about the essential idea of our 
factorization formula is in order. The basic idea of de- 
coupling the collinear and soft parts in the conventional 
factorization scheme is still respected here, but an im- 
portant step should be added to complete factorization. 
The separation of the collinear and soft parts is implicit 
in the conventional factorization scheme, however it is 
necessary to remove the soft limit in the collinear part 
properly at higher orders. It appears in loop calculations 
of the collinear part, and the contribution from the loop 
momentum in the soft (or usoft) limit should be removed 
since it is already included in the soft part. Our proce- 
dure of factorization exactly follows this idea. The soft 
contribution in the collinear part in loop calculations is 
extracted by the zero-bin subtraction, and there is no 
double counting of the soft degrees of freedom. 

As a result, the mixing of IR and UV divergences dis- 



appears, and the IR divergences are transformed to the 
UV divergences through the zero-bin subtraction as long 
as the zero-bin subtractions are performed properly. And 
the result for the kernel turns out to be IR finite. In this 
line of thought, the jet function in DIS takes a special 
place since it becomes IR finite after the zero-bin sub- 
traction. But considering the fact that the jet function is 
the matching coefficient, it is obvious that the jet func- 
tion is IR finite. 

We have shown that our factorization theorem factor- 
izes each part in a meaningful way by combining the soft 
part and the initial-state jet function. Since the initial- 
state jet function consists of the difference between the 
soft and usoft zero-bin contributions, the kernel W in 
its broad sense can be referred to as the soft part. The 
point is that the hard part and the kernel W are IR fi- 
nite, and can be computed in perturbation theory, and 
the evolution can be set up without concerns on how to 
treat IR divergences in each part. The detailed results 
will be presented in a forthcoming paper. This factoriza- 
tion method has been also applied to the processes with 
small transverse momentum . It strongly implies that 
our factorization idea can be widely used in various high- 
energy scattering processes. 
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